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Abstract 

A canonical variable coefficient nonlinear Schrodinger equation with a 
four dimensional symmetry group containing SL(2,M) group as a subgroup 
is considered. This typical invariance is then used to transform by a symmetry 
transformation a known solution that can be derived by truncating its Painleve 
expansion and study blow-ups of these solutions in the Lp-norm for p > 2, 
Loo-norm and in the sense of distributions. 
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1 Introduction 

It is well known that the linear heat and Schrodinger equations for u E M., ijj E C 

ut = Uxx, i'ipt + V'xx = 0, X eR, t > 

have isomorphic Lie symmetry groups. The symmetry group with the 
infinite- dimensional ideal reflecting linearity factored out can be written as a 
semidirect product of the three-dimensional Heisenberg group H and SL(2,M) 

G = H X SL(2,M). 
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On the other hand, among the modular class of nonlinear Schrodinger equations in 
one space dimension 

#t + V'xx = F(|7/'|)7/', (1.1) 

the only one preserving the symmetry group of the linear Schrodinger equation except 
with infinite-dimensional symmetry of its linear counterpart is the quintic Schrodinger 
equation 

#t + V^xx = A|7/^|V. (1.2) 

In this article we look at a variable coefficient extension of the one-dimensional 
cubic NLSE (nonlinear Schrodinger equation) 

i'il>t + il>xx + gix,t)\'ilj\'^'ip + h{x,t)ilj = 0, 

9 = 9i+i92, h = hi+ih2, gj,hj eR, j = 1,2, fi-i 7^ 0. 

Variable coefficient extensions of nonlinear evolution type equations tend to arise in 
cases when less idealized conditions such as inhomogeneities and variable topographies 
are assumed in their derivation. The reader is referred to [8] and the references 
therein for several physically motivated applications. Symmetry classification of 
(11. 4p was given in P]. We mention that the dimension of the maximal symmetry 
group is dimG = 5 and is achieved only when the coefficients are constant and 
additionally h = which is reduced to nothing more than the usual NLS equation. 
The symmetry group is isomorphic to the group of one-dimensional extended Galilei 
similitude algebra. We are particularly interested in the case when at least the 
SL(2,M) invariance is contained in the full symmetry group of (11.31) . This is usually 
referred to as the pseudo-conformal invariance in the context of qualitative analysis 
of PDEs. 

We emphasize that this invariance manifests itself as a subgroup of the full 
symmetry groups of the NLSE, Davey-Stewartson (DS) equations and their possible 
generalizations [5] and has been successfully applied to investigate blowup formation 
in these nonlinear evolution models [21 [TJ [H [3] . A study of self-similar solutions of 
the pseudo-conformally invariant nonlinear Schrodinger equation can be found in [6] . 
The point is that when the variable coefficients are allowed in these equations, this 
typical symmetry is mostly destroyed. Our intention here is to detect the subcases 
in which such a symmetry remains intact in variable coefficients variants of the NLS 
equations. This will make it possible to generate new nontrivial solutions from known 
ones. 

We quote the following result from |lj and note that throughout the paper, any 
solution will be understood as a pointwise solution in the classical sense. 

Proposition 1. Any equation of the form (ll.Sp containing SL{2,'R) symmetry group 
as a subgroup can be transformed by point transformations to the canonical form 

iiJt + i^xx + {e + i-f)-\ij\^ij + {hi+ih2)\iJ = 0, xeR\{0}, (1.4) 
where e = ±1, 7, hi and /i2 are arbitrary real constants. 
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Note that like all equations in the modular class (11. ip . equation fll.Sp and also 
fll.4p are always invariant under the constant change of phase of ip (gauge-invariance) 
while leaving the {x, t) coordinates unchanged. We represent the phase and modulus 
of Ip hj p, u writing ip = p{x, t) exp {iu^x, t)). 

Proposition 2. The symmetry algebra L of (11 ■4p is four- dimensional and spanned 
by the vector fields 



and W is the center element, namely commutes with all other elements. The Lie 
algebra L has the direct sum structure 



The elements T, D, C, W generate time translations, scaling, (pseudo)-conformal 
and gauge transformations, respectively. A significant consequence of these 
transformations is the group action on the solutions given by the following. 

Proposition 3. If4'o{x,t) is a solution of (11. 4p . then so is 



for ad — he = 1 . 

Proof. By exponentiating the infinitesimal generators T, D, C (i.e. solving Cauchy 
problems for these vector fields) and then composing the corresponding group 
transformations we find the above SL(2,M) action on the solutions. Note that the 
action corresponding to t generates Mobius transformations of t. 

Based on this result, our main purpose is to transform one known solution of the 
original equation to a more complicated one by the transformations of the SL(2,M) 
symmetry and then choose group parameters a, 6, c appropriately and pass to limit 
of the wave function ■?/' as t — )■ T~ for some finite time T in the Lp-norm and 
distributional sense as well. 

We now use a truncation approach to obtain a special exact solution of (II. 4p . We 
are going to find this special explicit solution by truncating its Painleve expansion at 
the first term. For convenience we write (ll.4p together with its complex conjugate as 
the system 




L = sl(2,M) ©M. 









0. 



(1.5) 
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Here u was employed instead of if) and v denotes its complex conjugate, but in this 
setting they are viewed as independent functions. We first show that (11. 5p does 
not pass the Painleve test for PDEs and then proceed to obtain an exact solution 
afterwards. 

A partial differential equation is said to have the Painleve property if all 
its solutions are single valued around any non-characteristic movable singularity 
manifold. If this singularity manifold is denoted by $(a;,t) = (actually a curve 
in this case), we shall look for solutions of the system (II. 5p in the form of a Laurent 
expansion and we expand 

oo oo 

u{x,t) = uj(x, t) (a;,t), v{x,t) = ^ Vj{x,t)<i>'^+\x,t), (1.6) 
i=o i=o 

where uo,vo ^ and Mj,fj,$(x, t) are analytic functions, a and /3 are negative 
integers to be determined from the leading order analysis, so as to ensure absence 
of essential singularities and branch points in all solutions. For the determination of 
leading orders a and /3, we substitute u ~ Mo$" and v ~ Vi^^^ in (II. 5p and see that 
by balancing the terms of smallest order 

a + /3 = -2 (1.7) 

and 

Mot^o = -«(«-!) ^-$' = -/3(/3-l)^-$' (1.8) 
e + «7 6 — ^7 

must hold. (II. 7p allows the negative integers a = — 1 and /3 = — 1 and with these 
leading orders, (II. 8p forces 7 = 0. After determination of the leading orders, we 
substitute (II. 6p into (II. 5p . Equating to zero the coefficient of j > 1, j G N, we 

arrive at a linear system 

«w (:;)-( ) 

from which the coefficients Wj, Vj can be found. Those values of indices j for which 
det(5(j) = are called resonances. In order that the expansion (II. 6p includes correct 
number of arbitrary functions as required by the Cauchy-Kovalewski theorem (where 
$(x,t) should be one of the arbitrary functions), some consistency conditions at 
resonances must be satisfied. In the general case (II. 3p . these constraints force the 
coefficients to be properly related (See [S] for details). For (II. 4p . this is not the case, 
namely it cannot pass the Painleve test. It thus fails to satisfy the necessary condition 
for the equation to have Painleve property. 

Notwithstanding this fact, application of the Painleve expansion to nonintegrable 
PDEs like (ll.4p (or more properly partial integrable) can allow particular explicit 
solutions to be obtained by truncating the expansion. This approach imposes 
constraints on the arbitrary functions and the function $ as a result of compatibility 
of an overdetermined PDE system. 
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Before truncating the series fll.6p at some index j = N, we first weaken the 
condition that a and /3 assumes only integer values. When we solve (11. 7p and (11.81) 
together, we find that for 7 7^ 



a 



-l-i6, f3 = -l + i6; 
and (II. 8p is equivalent to 
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-3e ± + 9 

2^ 



UqVq 



— 
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;i.io) 



;i.ii) 



Truncating the Painleve expansion at the first term (j = 0), we assume a solution of 
the form 



v{x,t) = vo{x,t)^{x,t) 



;i.i2) 



We substitute these in (II. 5p and set the coefficients of the terms $-3±i<5^ ^-2±i5^ <|)-i±*<5 
equal to zero. The condition at the order ^-^^^^ is equivalent to (11.110 and terms of 
order $-2±*<5^ ^-lif? disappear if $, uq, vq are chosen to be 



^{x,t) = 
Uo{x,t) = Ai 
vo{x,t) = A2 



X 



kit + ki 



2/3 



+ k2 



X 



1/6 



(M + fci)2/3 
^,1/6 

{kit + A;i)2/3 



exp 
exp 



A{kit + ki) 
kix"^ 
A{k4 + ki 



+ h 



;i.i3) 
;i.i4) 
;i.i5) 



for arbitrary real constants Ai, A2 and /ci, /C4 with the constraints hi = 5/36 and 
/12 = 0. We require to have v = u*, which implies Ai = A2. Let us rename this 
constant A. When we check the condition (II. lip , we find that = Since this 
square must be positive, it is necessary that 



6_ 
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-3e ± V872 + 9 



< 0. 



;i.i6) 



For both e = ±1, it is seen that the minus sign must be picked in the formula for S 
of (ll.lOp . Having found a consistent truncation, we can write the solution to (ll.4p as 



u{x, t) 



^^1/6 



exp 



A{kit + ki 



x^l^ + k2{kit + kifl^ 
Choosing ^4 = 0, we obtain the stationary solution 
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X 



2/3 



{kit + A:i)2/3 



U[X] 



x2/3 + kl'^k2 



exp 



5\n{x^'^ + kl''^k2) + h 



k2 \ +h 
(1.17) 

(1.18) 



where 



hi 



36^ 



;i.i9) 



and ki,k2, [k^ is relabelled) are arbitrary real constants. 
We summarize: 



Proposition 4. The following solves equation (11. 4p for arbitrary constants ki, k2, k^ 
and for the parameters hi = 5/36, h2 = 



u{x, t) = A' 



X 



1/6 



where A = {-^^/^ and 6 = z2lz^0fll. 



-dlnix"^/^ + kl^\2) + h 



;i.20) 



2 Transforming solutions by SL(2, M) group 

Blowup in the Lp, L^o norms and in the distributional sense 

Now we would like to illustrate how the SL(2,]R) group action can be useful in 
establishing blow-up profiles of initial value problems for variable coefficient NLS 
equations just as they were used for their constant coefficient counterparts. 

Let u{x) be the stationary solution to (II. 4p . defined by (I1.18I) - (I1.19I) . We set 
ipolx) := u{x) and use Proposition 3. By this proposition for arbitrary a, 6 G M 

V.(..t)^(a + «)-"'^exp(^)..„(-^) (2,1) 

is also a solution to equation (II. 4p . 

We can assume a > 0, 6 < and denote e := a + bt = b{t — ^) = b{t — T), where 
T = — I > 0. Hence, t — )■ <^ e — )■ 0+. By using this notation we can write solution 
(12. ip in the form 

V'.(x) =e-^^^exp (^j^o(-)- 

We are going to show that these solutions will blowup in the Lp-norm when p > 2, 
Loo-norm and in the sense of generalized functions, respectively. 



Note: In the following theorems, we do not impose an initial condition but rather 
we limit to the one as dictated by the solution (12. ip in the form 



^{x, 0) = 




Also, for a blow-up at some finite time we have to fix b figuring in (12.10 . 
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Lp-blow-up solutions. 

Theorem 1. For any T > there is a solution ilj{x,t) to equation (II ■4p such that 

lim ||^/'(x,t)|L = +CX) forallp>2, (2.2) 
i/p 



where \\'tp{x,t)\\p = {^J_^ \^\'p dx 

Proof. Let T > be a finite time. We can always arrange two numbers a > and 
6 < such that T = — |. Setting these numbers in (12. ip we get the function ilj{x,t) 
which will be instrumental in the proof. We rewrite this function, as given above, in 
the form 

^,(x)=.-V2exp(^)^o(f). 

Then we have 

\im Us{x)\\,= hm mx,t)\\p- (2.3) 
By the definition of ipoi^) we have 

V'o(-) = ^2 2 ^ ^ exp{iip{x,e)), 
^ e + 6 

/ 2 2 \ 2 

where ip{x, e) = — 5 ln( (|) ^ + k2 + k^) , C = kl /c2 and we choose C > 0. Thus, 

= — ^exp[i(— + <^{x,e))], 

X3 +e3C ^ \4:e J' 



and 



\ibAx^'^dx = 2 \ — ■ ■ — r-^dx. 

'^^ ^' Jo (x2/3+52/3C)P 



By (2.3) lim^^^r- ||-?/'(x, t) ||p = +oo if and only if: 

i) \4'e{x)\^ is finite for all p > 2 and e > 0, 

ii) lime_^o+ /o°° |V^e(x)p(ix = oo. 
Let us substitute x = ey. Then 



J — c 



\ij,{x)\Pdx = 2AP I 2P, 2 . ^,p ^dy = 7(^z^ / J^—^dy 



2AP 



r(p-2)/2 
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1 ^ 

where j-^ — yd-y is convergent by continuity of the function. On the other hand, 



s. P 
ye ye 



< 



[yi+C) yif 1/2 



p p 



Consequently, ^ 2^ y dy < -^dy imphes that the integral j-r^-^dy is 

convergent for all p > 2. Hence, \ipe{x)\^dx is convergent for p > 2. Finally, by 
taking the limit we find 

lim [ \ilj,{x)\Pdx = lim [ 2^\.p dy = 00. 

Blow-up solutions in L^o- norms. In the following theorem we prove that the 
above defined solutions ipelx) will blowup in Lqo- norm too. 

Theorem 2. For any T > there is a solution ip{x,t) to equation (11. 4p such that 

lim ||V^(x,t)||oo = 00, (2.5) 

where \\'ib{xA)\\c^ = esssup^^rn^^ l-il^ixAl, t <T. 



ess sup2.g[o,oo) 

Proof. Let again ip{x,t) be the function defined by (12. ip and ipe{x) be the above 
defined function. By the construction of ipelx) 

lim \\ipe{x)\\oo = lim \\'ip{x,t)\\ 00 

and therefore we are done if we can show that lime^o+ 1 1 V^e (3^)11 00 = 00. We have 
|'V^e(3^)| = — T^^^T— • Let A = C = 1. \ipi;{x) \ is an even function so that we can restrict 
ourselves to the interval [0, 00). Since |'?/'e(a;)| is continuous on [0, 00), |V^e(0)| = and 
lim;j._j.oo |V^e(2;)| = 0, then there exists Xq G (0, 00) such that 



\\ipe{x)\\oo = ess sup \ilJe{x\= max \ipe{x)\ 



Xq 



xe[o,c») ^6[o,oo) x^+elC 



A simple computation yields xq = Therefore, 



Consequently, 



lim ||V'e(x)||oo = 00. 

£-5-0+ 
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5-Blowup solutions in the sense of generalized functions. Let D = C^(0, oo) 
be the space of infinitely differentiable functions with compact support in (0, oo). The 
dual of D is called the space of generalized functions and is denoted by D'. 

Definition 1. Let f, fn G D' . We say that the sequence fn converges to f if and only 

^/ {fn,f) {J\f) for all if E D , where {f,ip) denotes the value of the functional f 
at ip. 

Now we present 5-blowup solutions in the sense of generalized functions. 
Theorem 3. For all p > 2 

e(p-^)/^\^^(^x)\p ^ KS{x) ase^O+, 

oo I I K 

where K = -r-^ — srdy, C > and 5(x) denotes the Dirac distribution at the 
J-^ [yi+c) ^' ^ ' 

origin. 

Proof. By Definition [T] we have to show that 

/oo 
e^P-'^)/^\^Jx)Y'^{x)dx = Kip(0) for all cp e D. 
-oo 

Evidently, 



e 



(^-^)/^/ \Mx)\Mx)dx = e^^-'y' I ^j-^^^^^^{x)dx. 



By setting the substitution x = ey we obtain 

r mx)\Mx) dx = e^^-^y^ r i^Zl''^Z^p ^iey)edy 



e 



^^-^^'^ I , f}'^:" ,Mey) dy - r T^T^^iey) dy. 



(2.6) 



Let feiy) := ,^^'^^'^^\p 'p{ey). The sequence feiy) satisfies the following two conditions: 

i) \fe{y)\ < Cs-^^ e Li(-oo, oo) if p > 2, 

ii) lim,^o+/.(y) = -^^^^'^(0). 

Then by Lebesgue's dominated convergence theorem we obtain 

/oo /.oo y^PlylP/S 

fMdy = m / ( 2,1] = ^^^0)- 
-oo J-oo [y^'-^ + C) 

By Definition 1, this implies that 
as e — i- 0+. 

Remark: The above argument indicates that the singular behavior of the solution 
is even much worse than in the usual distributional sense as was done in Ref. [7]. 
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